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Abstract 

This paper investigates the performance of the point-to-point multiple-input-multiple-output (MIMO) 
systems in the presence of a large but hnite numbers of antennas at the transmitters and/or receivers. 
Considering the cases with and without hybrid automatic repeat request (HARQ) feedback, we determine 
the minimum numbers of the transmit/receive antennas which are required to satisfy different outage 
probability constraints. Our results are obtained for different fading conditions and the effect of the power 
amplihers efficiency on the performance of the MIMO-HARQ systems is analyzed. Moreover, we derive 
closed-form expressions for the asymptotic performance of the MIMO-HARQ systems when the number 
of antennas increases. Our analytical and numerical results show that different outage requirements can 
be satished with relatively few transmit/receive antennas. 

I. Introduction 

The next generation of wireless networks must provide data streams for everyone everywhere 
at any time. Particularly, the data rates should be orders of magnitude higher than those in 
the current systems; a demand that creates serious power concerns because the data rate scales 
with power monotonically. The problem becomes even more important when we remember 
that currently the wireless network contributes ~ 2% of global CO 2 emissions and its energy 
consumption is expected to increase 16 — 20% every year lH]. 

To address the demands, the main strategy persuaded in the last few years is the network 
densification [l2l . One of the promising techniques to densify the network is to use many antennas 
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at the transmit and/or receive terminals. This approach is referred to as massive or large multiple- 
input-multiple-output (MIMO) in the literature. 

In general, the more antennas the transmitter and/or the receiver are equipped with, the better 
the data rate/link reliability. Particularly, the capacity increases and the required uplink/downlink 
transmit power decreases with the number of antennas. Thus, the trend is towards asymptotically 
high number of antennas. This is specially because millimeter wave communication Q, Hi, 
which are indeed expected to be implemented in the next generation of wireless networks, makes 
it possible to assemble many antennas at the transmit/receive terminals. However, large MIMO 
implies challenges such as hardware impairments and signal processing complexity which may 
limit the number of antennas in practice. Also, one of the main bottlenecks of large MIMO is 
the channel state information (CSI) acquisition, specifically at the transmitter. Therefore, it is 
interesting to use efficient feedback schemes such as hybrid automatic repeat request (HARQ) 
whose feedback overhead does not scale with the number of antennas. 

The performance of HARQ protocols in single-input-single-output (SISO) and MIMO systems 
is studied in, e.g., O-lIHl and Il9ll- [l20ll . respectively. MEMO transmission with many antennas is 
advocated in [|2T1, ll^ where the time-division duplex (TDD)-based training is utilized for CSI 
feedbaclli. Also, [|^-[|29l introduce TDD-based schemes for large systems. In the meantime, 
frequency division duplex (FDD)-based massive MIMO has recently attracted attentions and 
low-overhead CSI acquisition methods were proposed Considering imperfect CSI, 

[1^ derives lower bounds for the uplink achievable rate of the MIMO setups with large but 
finite number of antennas. Finally, lf34]l (resp. [|35l) studies the zero-forcing based TDD (resp. 
TDD/FDD) systems under the assumption that the number of transmit antennas and the single¬ 
antenna users are asymptotically large while their ratio remains bounded (For detailed review of 
the literature on massive MIMO, see lf^ - [[38ll l. 

To summarize, a large part of the literature on the point-to-point and multi-user large MIMO 
is based on the assumption of asymptotically many antennas. Then, a natural question is how 
many transmit/receive antennas do we require in practice to satisfy different quality-of-service 
requirements. The interesting answer this paper establishes is relatively few, for a large range of 

'The results of CQ-OEI are mostly on multi-user MIMO networks, as opposed to our work on point-to-point systems. 
However, because many of the analytical results in EQ-OHl are applicable in point-to-point systems as well, these works are 
cited. 
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outage probabilities. 

Here, we study the outage-limited performanee of point-to-point MIMO systems in the eases 
with large but finite number of antennas. We derive elosed-form expressions for the required 
number of transmit and/or reeeive antennas satisfying various outage probability requirements 
(Theorem 1). The results are obtained for different fading eonditions and in the cases with or 
without HARQ. Furthermore, we analyze the effect of the power amplifiers (PAs) efficiency and 
the antennas spatial correlation on the system performance (Sections IV.A and V.B, respectively) 
and study the outage probability in the cases with adaptive power allocation between the HARQ 
retransmissions (Section IV.B). Finally, we study the asymptotic performance of MIMO systems. 
Particularly, denoting the outage probability and the number of transmit and receive antennas 
by Pr(Outage), At and At, respectively, we derive closed-form expressions for the normalized 
outage factor which is defined as F = — i°g(Pd^^tage)) number of transmit and/or receive 

antennas increases (Theorem 2). 

As opposed to ifSll- EOll . we consider large MIMO setups and determine the required number 
of antennas in outage-limited conditions. Also, the paper is different from [l2T]| - [[^ because we 
study the outage-limited scenarios in point-to-point systems, implement HARQ and the number 
of antennas is considered to be finite. The differences in the problem formulation and the channel 
model makes the problem solved in this paper completely different from the ones in [[5l- l[38ll . 
leading to different analytical/numerical results, as well as to different conclusions. Finally, our 
discussions on the asymptotic outage performance of the MIMO setups and the effect of PAs 
on the performance of MIMO-HARQ schemes have not been presented before. 

Our analytical and numerical results indicate that: 

• Different quality-of-service requirements can be satisfied with relatively few transmit/receive 
antennas. For instance, consider a SIMO (S: single) setup without HARQ and transmission 
signal-to-noise ratio (SNR) 5 dB. Then, with the data rate of 3 nats-per-channel use (npcu), 
the outage probabilities 10“^, 10“^ and 10“^ are guaranteed with 16, 18 and 20 receive 
antennas, respectively (Fig. 5a). Also, the implementation of HARQ reduces the required 
number of antennas significantly (Fig. 3). 

• Considering the moderate/high SNRs, the required number of transmit (resp. receive) an¬ 
tennas scales with (Q“^(0))^, and linearly, if the number of receive (resp. 

transmit) antennas is fixed. Here, M is the maximum number of HARQ retransmission 
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rounds, T is the number of ehannel realizations experienced in each round, (p denotes the 
SNR, 9 is the outage probability constraint and represents the inverse Gaussian Q- 

function. These scaling laws are changed drastically, if the numbers of the transmit and 
receive antennas are adapted simultaneously (see Theorem 1 and its following discussions 
for details). 

• For different fading conditions, the normalized outage factor T = converges 

to constant values, unless the number of receive antennas grows large while the number of 
transmit antennas is fixed (see Theorem 2). Also, for every given number of transmit/receive 
antennas, the normalized outage factor increases with MT linearly. 

• There are mappings between the performance of MIMO-HARQ systems in quasi-static, 
slow- and fast-fading conditions, in the sense that with proper scaling of the channel 
parameters the same outage probability is achieved in these conditions. This point provides 
an appropriate connection between the papers considering one of these fading models. 

• Adaptive power allocation between the HARQ retransmissions leads to marginal antenna 
requirement reduction, while the system performance is remarkably affected by the ineffi¬ 
ciency of the PAs. Finally, the spatial correlation between the antennas increases the required 
number of antennas while, for a large range of correlation conditions, the same scaling rules 
hold for the uncorrelated and correlated fading scenarios. 

II. System Model 

Consider a point-to-point MIMO setup with Nt transmit antennas and A), receive antennas. 
We study the block-fading conditions where the channel coefficients remain constant during 
the channel coherence time and then change to other values based on their probability density 
function (PDF). In this way, the received signal is given by 

Y = BX + Z,ZeC^^^\ (1) 

where H G is the fading matrix, X G is the transmitted signal and Z G denotes 

the independent and identically distributed (IID) complex Gaussian noise matrix. The results are 
mainly given for IID Rayleigh-fading channels where each element of the channel matrix H 
follows a complex Gaussian distribution CA/'(0,1) (To analyze the effect of the antennas spatial 
correlation, see Fig. 7 and Section V.B). The channel coefficients are assumed to be known at 
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the receiver which is an acceptable assumption in block-fading channels If5l- ll^ . On the other 
hand, there is no CSI available at the transmitter except the HARQ feedback bits. The feedback 
channel is supposed to be delay- and error-free. 

As the most promising HARQ approach leading to highest throughput/lowest outage probabil¬ 
ity [El, lEl, im, l[T4l . we consider the incremental redundancy (INR) HARQ with a maximum of 
M retransmissions, i.e., the message is retransmitted a maximum of M times. Note that setting 
M = 1 represents the cases without HARQ, i.e., open-loop communication. Also, a packet 
is defined as the transmission of a codeword along with all its possible retransmissions. We 
investigate the system performance for three different fading conditions: 

• Fast-fading. Here, it is assumed that a finite number of channel realizations are experienced 
within each HARQ retransmission round. 

• Slow-fading. In this model, the channel is supposed to change between two successive 
retransmission rounds, while it is fixed for the duration of each retransmission. 

• Quasi-static. The channel is assumed to remain fixed within a packet period. 

Fast-fading is an appropriate model for fast-moving users or users with long codewords compared 
to the channel coherence time Id, lEl. On the other hand, slow-fading can properly model 
the cases with users of moderate speeds or frequency-hopping schemes llT0l - |[T^ . ifTSl - lfTTl . 
Finally, the quasi-static represents the scenarios with slow-moving or stationary users, e.g., [El, 

m, d-d. 


HI. Problem Formulation 

Considering the INR HARQ with a maximum of M retransmissions, Q information nats are 
encoded into a parent codeword of length ML channel uses and the codeword is divided into 
M sub-codewords of length L. In each retransmission round, the transmitter sends a new sub¬ 
codeword and the receiver combines all signals received up to the end of that round. Thus, the 
equivalent rate at the end of round m is — npcu where R denotes the initial transmission 

rate. The retransmissions continue until the message is correctly decoded by the receiver or the 
maximum permitted retransmission round is reached. 

Let us denote the determinant and the Hermitian of the matrix X by |X| and X*’, respec¬ 
tively. Assuming fast-fading conditions with T independent fading realizations H((m — 1)T + 
1),... ,H(mT) in the mth round and an isotropic Gaussian input distribution over all transmit 
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antennas, the results of, e.g., [|40l Chapter 15], IHTI Chapter 7], can be used to find the outage 
probability of the INR-based MIMO-HARQ scheme as 



M 


nT 






R 


< — • ( 2 ) 
- M ' ^ ^ 


Pr(Outage)^"‘"§ = Pr 

I ivi ; 

n=l t={n—l)T+l 

Here, 0 is the total transmission power and is the transmission power per transmit antenna 
(in dB, we have 10 log^g ^ which, because the noise variance is set to 1, represents the SNR as 
well). Also, Iat, represents the x identity matrix. 

Considering T = 1 in dH), the outage probability is rephrased as 


M 


Pr(Outage)^“‘^‘"s = M 


n=l 




< 


M 


( 3 ) 


in a slow-fading channel. Also, setting H(f) = H, Vf = 1,..., MT, the outage probability in a 
quasi-static fading channel is given by 




R 

< — 
- M 


( 4 ) 


Pr(Outage)Q“"^'-^“‘‘^ = Pr (^log 

Using ©-(HI) for given initial transmission rate and SNR, the problem formulation of the paper 
can be expressed as 

{iVt, iVr} = argmin{Pr(Outage) < 6}. (5) 

Nt,Nr 

Here, 9 denotes an outage probability constraint and Ni, N^- are the minimum numbers of 
transmit/receive antennas that are required to satisfy the outage probability constraint. In the 
following, we study ® in four distinct cases: 

• Case 1: is large but At is given. 

• Case 2: Ar is given but At is large. 

• Case 3: Both At and Ar are large and the transmission SNR is low. 

• Case 4: Both At and Ar are large and the transmission SNR is high. 

It is worth noting that the three first cases are commonly of interest in large MIMO systems. 
However, for the completeness of the discussions, we consider Case 4 as well. Moreover, in 
harmony with the literature na, eaB we analyze Cases 3-4 under the assumption 


Ar 


= A, 


( 6 ) 


^In 1341 . 1351 . which study multi-user MIMO setups, Nt and Ni are supposed to follow ® while, as opposed to our work, 
they are considered to be asymptotically large. 


May 1, 2015 


DRAFT 










7 


with K being a constant. However, as seen in the following, it is straightforward to extend the 
results of the paper to the cases with other relations between the numbers of antennas. 


IV. Performance Analysis 


To solve (|5]), let us first introduce Lemma 1. The lemma is of interest because it represents 
the outage probability as a function of the number of antennas, and simplifies the performance 
analysis remarkably. 

Lemma 1: Considering Cases 1-4, the outage probability of the INR-based MIMO-HARQ 
system is given by 

Vr(Outage)P^“"g = Q ^ (i) 

< Pr(Outage)= Q , (ii) (7) 

Pr(Outage)Q“"^‘-^‘^‘‘‘^ = Q (^) , (iii) 


where Q{x) = e~^du is the Gaussian Q-function and for different cases fi and a are 

given in dS]). 

Proof: The proof is based on dH)-® and |l42l Theorems 1-3], where considering Cases 
1-4 the random variable Z(t) = log \1n, + converges in distribution to a Gaussian 

random variable Y ~ AA(/i, cr^) which, depending on the numbers of antennas, has the following 
characteristics 


= < 




At 

Nr J ’ 


A = iV™.log (a) + At..,,. (eS 


if Case 1 
if Case 2 
if Case 3 
if Case 4 


Njnax—Nmin 1 


7 


E Nmin-l 
i=l 


Nrr 


—, 7 = 0.5772... 

X-* ’ ' 


a 


i=l 


N —1 

-‘’'min 


+ N„. 


Nmsx = max(At, N,), V^in = min(At, N,). 


E -/Vma 
2 = 1 


-^max ~ 1 1 


( 8 ) 


In this way, from dH) and for different cases, the outage probability in fast-fading condition is 
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\ 1 
' t=l 

where, because Z is the average of MT independent Gaussian random variables Y ~ AA(/i, a^), 
we have Z ~ A/”(/i, Consequently, using the cumulative distribution function (CDF) of 

the Gaussian random variables, the outage probability of the fast-fading condition is given by 
dlli). The same arguments can be applied to derive ([Tlii-iii) in the slow-fading and quasi-static 
conditions. ■ 

Lemma 1 leads to the following corollaries: 

1) For Cases 1-4, using INR MIMO-HARQ in the quasi-static, slow- and fast-fading channels 
leads to scaling the variance of the equivalent random variable by 1, M and MT, respec¬ 
tively. That is, using HARQ, there exists mappings between the quasi-static, the slow- and 
the fast-fading conditions in the sense that with proper scaling of a in (|7]) they lead to the 
same outage probability. 

2) With asymptotically large numbers of transmit and/or receive antennas, the optimal data 
rate which leads to zero outage probability and maximum throughput is given by i? = 
11 — io,u 0, with fi derived in dS]); Interestingly, the result is independent of the fading 
condition. Also, with asymptotically high number of antennas and R = fi — u,io^0, no 
HARQ is needed because the message is decoded in the first round (with probability 1). 

3) Finally, using ([2]), we can map the MIMO-HARQ system into an equivalent SISO-HARQ 
setup whose fading realizations follow AA(/U., with /x and a given in ([8]) for different 
cases. 

Using Lemma 1, the minimum numbers of antennas satisfying different outage probability 
constraints are determined as stated in Theorem 1. 

Theorem 1: The minimum numbers of the transmit and/or receive antennas in an INR-based 
MIMO-HARQ system that satisfy the outage probability constraint Pr(Outage) < 9 are given 


Pr(Outage)“‘^‘"§ = Pr ( Z < ^ 
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(pVTTrQ-He) 


/MTil+<t>){Nrlog(l+<j>)- 


N, = N,N, = KN,N = ^ + 


_ K , Q~He) 


M(j> ^ s/MTK’’ 


N, = N,Nt = KN 
N ~ 

iV~ 




log((/))- 7 -l + (/4'-l)log(-j^) ’ 


if Case 1 

if Case 2 

if Case 3 
if Case 4 

K > 1 

iC < 1, 


(10) 


i^(log( 0 )- 7 -l-log(A')+(i^)log(l-i^)) 

if the channel is fast-fading. Here, Q~^{x) and W{x) denote the inverse Q-function and the 
Lambert W function, respectively. For the slow-fading and quasi-static conditions, the minimum 
numbers of the antennas are obtained by (fT^ where the term -—— 

Q~^{9), respectively. 


/MT 


is replaced by and 


Proof: Considering Lemma 1 and a fast-fading condition, ([5]) is rephrased as 

A(^)l 


{Nt, = arg 

Nt,Nr 


F M 

a 


which for different cases leads to 
Case 1: 


N, = arg <( iVt log ( 1 -f 

Nr 


N4 


_R 

M 


VWT J ’ 


Q-\o) [K 


( 11 ) 


(“) , 

~ — arg < log(w) = 


Nt J M y/Mf V N, 
R ^ Q-\o)V^ \ 


Nr = 


MN, y/MTN,^ J 


AMTN,W^ 


( 12 ) 


Case 2: 


Nt = arg <j Nr log(l + f) - — = 


R Q-\e)fy/K \ 


^ Nt = 

Case 3: Nr = N, Nt = KN where 


M y^MT{l + f)y/Nt! 

fVKQ-Ko) 


VMT{1 + <P) (iV,log(l + 0)-f) 


N = sugiNf - — = 

n \ M y/MTK J 


N = ^+ 


M0 ^MTK' 


(13) 


(14) 
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In (fT^ . (a) is obtained by using the approximation log(l + u) ~ log(M) for large u’s and variable 
transform u = Also, W{x) denotes the Lambert W funetion defined as ye^ = x ^ y = 
W{x) |l43l. Note that the Lambert W funetion has an effieient implementation in MATLAB and 
MATHEMATICA. 


For Case 4, we eonsider two seenarios and use the following approximations. 
Case 4 with K > 1: Then, N^. = N,Nt = KN and 


N = arg N log 

N 


KN 


{K-l)N 


N-1 


+ E T-T +E 


R 


i=l 


2=1 


KN-i M 


N-l 

E 


Q-\o) _ 

y/MT\\i^, {{K-l)N + if ' \6 


KN-l 




2 = 1 


W . 

~ arg < N log 


N 


KN 
Q-\0) 


+ N (log((A: - 1)A^) - 7 ) + 2 - AT + a: log 


KN-l 
{K-l)N + l 


M 


(AT - 1)A^(2 - N) 


./MT \ {KN -N + 1){KN - 1) 


log 


(c) 


N 


N = 


~ arg < N log(0) — 7 — 1 + (K — 1) log 


K 


K-1 


KN-l 
{K-1)N + 1 

^ Q-\9) 
M ./WT 


KN-l 




2=1 


'log 


K 


K-1 


^ Q-^e) 


M y/MT 


ilh) 


log((/)) - 7 - 1 + (a: - 1) log {-^) ■ 


(15) 


Here, {b) is obtained by implementing the Riemann integral approximation /(*) — 

in the first three summation terms. Then, (c) follows from some manipulations, the faet that N 

is assumed large, and N{'^ — ^ R large A^’s. 
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Case 4 with K <1\ Then, Nj. = N,Ni = KN and 


N = arg KN log 

N 


KN 


(l-K)N 


KN-1 


+ ™ E 7-^ + E 


R 


2 = 1 


2 = 1 


N-i M 


Q-\0) 


KN-l 

E 


2 1 


VMT\ ^ {{l-K)N + iy 




2=1 


(fe) 

~ arg <( KA^log 

N 


KN 


+ KN (log((l - K)N) -j) + 2- KN + N\og 


N-1 


{1-K)N + 1 


R 

M 


Q-HO) 


(1 - K)N{2 - NK) 


y/WT\ {N-NK + 1){N-1) 


log 


A^- 1 


AT-l 


T-Ei 


I. 


N 




~ arg NK ( log((/)) — 7 — 1 — log(ii') + 


K-1 

K 


{1-K)N + 1^ 

log(l - ^)) - I = ^V-log(l -/0 


N = 


i +wv'-Mi-a') 


(16) 


K (log(^) - 7 - 1 - log(/f) + (V) log(l - K))' 
where ( 6 ) and (c) are obtained with the same procedure as in I117II1 . Finally, note that with 
slow-fading and quasi-statie fading ehannel (fTTT) is rephrased as 


and 


{iVt, = arg 

Nt,N, 


{iVt, = arg 

Aft.Afr 


F - f _ Q-HO) 


a 


a — — 
a 


\/M 


= Q-m , 


(18) 


(19) 


respeetively. Therefore, as stated in the theorem, with slow-fading and quasi-statie ehannel the 
required numbers of the transmit and/or reeeive antennas are determined by (fTO while the term 




is replaeed by 


Q-he) 


and Q ^( 6 ^), respeetively. 


Aeeording to Theorem 1, the following eonelusions ean be drawn: 


^We can follow the same procedure as in GUl-Glll to write 

^ |^(iog(i/’) ~ ~ ~ ^ -1) +1|. (17) 

in the cases with K = 1 which can be solved numerically via, e.g., “’fsolve” function of MATLAB or by different approximation 
schemes. However, for simplicity and because it is a special condition, we do not consider if = 1 as a separate case. 
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1) Using the tight approximation W(e°'~^^) ~ a: + a — log(a + x) in (fTOl) . the required number 
of receive antennas in Case 1 is rephrased as 


Nr ~ 




4MTNr (log 


R 


2MNt 


- log log 


Q-He)V^ 


R 


2MNi 


o, (20) 

MTiVt (log(0))' 

where the last approximation holds for moderate/high SNRs. Thus, at moderate/high SNR 
regimes, the required number of receive antennas increases with {Q~^{9)Y linearly. On the 
other hand, the required number of receive antennas is inversely proportional to the number 
of experienced fading realizations MT, the number of transmit antennas A^t, and (log(0))^. 
Interestingly, we can use (fTOlCase 2) to show that at high SNRs the same scaling laws hold 
for Cases 1 and 2. That is, in Case 2, the required number of transmit antennas decreases 
(resp. increases) with MT, A), and (log(^))^ (resp. {Q~^{9)y) linearly. 

2) The same scaling laws are valid in Cases 3 and 4, i.e., when the numbers of transmit and 
receive antennas increase simultaneously. For instance, the required number of antennas 
increases with Q~^{9) and the code rate R semi-linearlj^ (see (fTOlCases 3-4)). At hard 
outage probability constraints, i.e., small values of 9, the required number of antennas 
decreases with the number of retransmissions according to On the other hand, the 
number of antennas decreases with M linearly when the outage constraint is relaxed, i.e., 
9 increases. The only difference between Cases 3 and 4 is that in Case 3 (resp. Case 4) the 
number of antennas decreases with (f) (resp. log(()))) linearly. 

3) It has been previously proved that at low SNRs the same performance is achieved by the 
MIMO systems using INR and repetition time diversity (RTD) HARQ [[T4l Section V.B]. 
Thus, although the paper concentrates on the INR HARQ, the same number of antennas 
are required in the MIMO-RTD setups, as long as the SNR is low. 

As the number of antennas increases, the CDF of the accumulated mutual information, e.g., 
A?T ^n=i Yjt={n-1)T+1 + ;|[H(t)H(f)''| in fast-fading conditions, tends towards the step 

function. Therefore, depending on the SNR and the initial rate, the outage probability rapidly 
converges to either zero or one as the number of transmit and/or receiver antennas increases. To 


''The variable y is semilinear with a; if 1 / = a -|- foa; for given constants a and b. 
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further elaborate on this point and investigate the effeet of the number of antennas, we define 
the normalized outage faetor as 

log(Pr(Outage)) 


r = 


iVtiVr 


( 21 ) 


Intuitively, (1^ gives the negative of the slope of the outage probability curve plotted versus 
the product of the numbers of transmit/receive antennas. Also, (1211) follows the same concept 
as in the diversity gain D = — lim<^^oo [IH Eq. 14] which is an efficient metric for 

the asymptotic analysis of the MIMO setups. Theorem 2 studies the normalized outage factor 
in more details. 

Theorem 2: For Cases 1-4, different fading conditions and appropriate initial rates/SNR, the 
normalized outage factor is approximated by (l24l) . 

Proof: With given initial transmission rate R and SNR we use the approximation 


Qi^) - 0, 

for large x’s and (|7|) to rewrite the normalized outage factor (1211) as 


( 22 ) 


r = 




2N,N,a^ ’ 

ME, Fast-fading 
c = { M, Slow-fading 
1, Quasi-static. 

Then, from dSj), the normalized outage factor in different cases is found as 


(23) 


c 1 


2NiNr 

2 

c 

(Wlog(l-F0)-^) N,{l+^f 

2NtNr 


c 

c 

2KN^ 

^ 2 

K 

—a 



Nr<l> 

Nt 


R 

MNt 


2 


(log(l + f) 


R 

MNr 


if Case 1 
if Case 2 
if Case 3 
if Case 4 


r = 


a = 


(log(</,)+(K-l)log(^)-7-l) ^ 


(24) 


jr^(iog(0)-i-.s^ iog(i-rr)-iog(jr)-7-i) . 


-iog(i-ir) 

ME, Fast-fading 
c = < M, Slow-fading 
1, Quasi-static, 


if K < 1, 
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where a is found by following the same approaeh as in (fT5l)-(fT^. Finally, note that to use (l2^ 
the initial rate and the SNR should be sueh that > f for the eonsidered number of antennas. 
Otherwise, the outage probability eonverges to 1 and F —)■ 0. ■ 

Interestingly, the theorem indieates that: 

1) The normalized outage faetor beeomes eonstant in all eases, exeept Case 1 with a given 
(resp. large) number of transmit (resp. reeeive) antennas. Intuitively, this is beeause in all 
eases (exeept Case 1) the power per transmit antenna deereases by increasing the number 
of transmit antennas. Therefore, there is a tradeoff between increasing the diversity and 
reducing the power per antenna and, as a result, the normalized outage factor converges to 
the values given in (l2^ . In Case I, however, the message decoding probability is always 
increased by increasing the number of receive antennas and, as seen in Theorem 2, the 
normalized outage factor increases with monotonically, as long as fi > ^. 

2) In Case 3, the normalized outage factor becomes independent of the transmission SNR as 
long as /i > ;^. In Cases I, 2 and 4, on the other hand, the normalized outage factor scales 
with the SNR according to (log(0))^, if the SNR is high. 

3) In all cases, the normalized outage factor scales with the number of experienced fading 

realizations during the packet transmission, i.e., MT, linearly. Note that the same conclusion 
has been previously derived for the diversity gain D = — lim^^oo ifTOl . ifT^ . 

4) In cases 3-4, the normalized outage factor does not depend on the initial transmission rate. 
Moreover, in Case 3 the normalized outage factor is independent of the ratio between the 
number of transmit and receive antennas. 

A. On the Effect of Power Amplifiers 

As the number of the transmit antennas increases, it is important to take the efficiency of 
radio-frequency PAs into account Il3^ - ll38ll . For this reason, we use Lemma 1 to investigate the 
system performance in the cases with non-ideal PAs as follows. 

It has been previously shown that the PA efficiency can be written as Il44]l . Il45ll . Il46l Eq. (3)] 
and Eq. (3)] 

0cons \ 0max J 




(25) 
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Here, (f), and are the output, the maximum output, and the eonsumed power of the PA, 
respeetively, e G [0,1] denotes the maximum power effieieney aehieved at 0 = and 'd is 
a parameter that, depending on the PA elasses, varies between [0,1]. In this way, and beeause 
the INR-based MIMO-HARQ setup ean be mapped into an equivalent SISO-HARQ system (see 
Lemma 1 and its following diseussions), the equivalent mean and varianees ([8]) are rephrased as 


' (iV,log(l + f-y^).f). 

Nr log (l + ^ 


\ 


Ml 1+ 
2 


1_,5/ (,/,max)i 5 




j\T M f 

-'''r Y ; JV, Y 


Case 1 


, Case 2 


Case 3 


Case 4 


A = Nrrrrr log (EtT'"- i “ oj + EST , 

~2 —1 i \ AT f V^Mmax —1 1 \ 

^ - 2^i=l (7V™ax-M„,„+i)2 g -2 J , 

Amax = max(At, AJ, A^min = min(At, A^), 7 = 0.5772 .... 


(26) 


in the eases with non-ideal PAs. This is the only modifieation required for the non-ideal PA 
seenario and the rest of the analysis remains the same as before. 

B. On the Effect of Power Allocation 

Throughout the paper, we studied the system performanee assuming a peak power eonstraint 
at the transmitter. However, the system performanee is improved if the transmission powers are 
updated in the HARQ retransmission rounds. 

Let the transmission power in the mth round be fm- Then, the outage probability in the 
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fast-fading conditioi]^, i.e., Q, 
Pr(Outage)“‘*“§ = Pr 


is rephrased as 

M mT 


Bpr 


/^(m) 



m=l t=(m-l)T+l 
M 

Y.z„ 

m=l 




< 


R 

M 


1 


m 





n=l 




(27) 


n=l 


where fin and cr„ are obtained by replaeing into ([8]). Here, (d) is obtained by = 
log + ^H(f)H(f)^| ~ Af{Tfim,Ta'^). Also, (e) is based on the faet that the 
sum of independent Gaussian random variables is a Gaussian random variable with the mean 
and varianee equal to the sum of the variables means and varianees, respeetively. 

If the message is eorreetly decoded in the mth round, the total transmission energy and the 
total number of channel uses are ^(m) = and /(m) = rnL, respectively. Also, the total 

transmission energy and the number of channel uses are = -^ J2n=i 4 >m and I^m) = ML if an 
outage occurs, where all possible retransmission rounds are used. Thus, we can follow the same 
procedure as in jSll, (ill, (141 to find the average power, defined as the expected transmission 
energy over the expected number of channel uses, as 


$ = 


+ Em=l Pr EEi J2t={n-l)T+l 

I* + 

< 

r ) 

m j 

1 + EL.‘ (ife EEi EEL-dt+i i°g 


VI 





(m) 


1 + E™=i Q 


M(m) 


Z{m) 


(28) 


In this way, with a power constraint $ < 0, the problem formulation (5]) is rephrased as 

{iVt,7Vr} = argmin{Pr(Outage) < 9} 

Nt,N, 

s.t.i> < (j), 

which, using (8]) and (ITT]) , can be solved numerically or analytically. 


(29) 


^For simplicity, the results of this part are given mainly for the fast-fading condition. It is straightforward to extend the results 
to the cases with other fading models. 
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V. Simulation Results and Discussions 

In this section, we verify the accuracy of the derived results, and present the simulation results 
in spatially independent and correlated fading conditions as follows. 

A. Performance Analysis in Spatially-independent Fading Conditions 

In Figs. 1-4, we verify the accuracy of the results in Theorem 1 and derive the required 
number of transmit/receive antennas in the outage-limited conditions. Setting M = 2, = 1 

(Case 1), and the outage probability constraints Pr (Outage) < 9 (with 9 = 10“^, 10“^), Fig. 1 
shows the required number of receive antennas versus the initial transmission rate R. The results 
of the figure are obtained for the slow-fading conditions and different transmission SNRs. Then, 
considering A^r = 1 or 2, Fig. 2 demonstrates the required number of transmit antennas in Case 
2 with large and given N^. Here, we consider quasi-static, slow- and fast-fading conditions 
with 9 = 10“'^, T = 2, M = 2, 0 = 15 dB. In Fig. 3, we verify the effect of HARQ on the 
system performance. Here, assuming Case 1 (large and iV) = 1,5), the required number of 
antennas is derived in the scenarios with (M = 2) and without (M = 1) HARQ. The results of 
the figure are given for quasi-static channels, ^ = 5 dB and 9 = 10“'^. 

Figure 4 studies the required number of antennas in Cases 3 and 4 with low and high SNRs, 
respectively, large number of transmit and receive antennas, and ^ = K. Also, the figure 
demonstrates the analytical results of Theorem 1 when the approximation steps (b) — (c) of 
(fT5])-(fT^ are not implemented, i.e., (fTTl) is solved numerically via dU). Here, we consider the 
quasi-static conditions, M = 1, and 9 = 10“^. Note that, to have the simulation results of Case 
4 in reasonable running time, we have stopped the simulations at moderate initial transmission 
rates. For this reason, the simulation results of Case 4, i.e., the red solid-line curves of Case 4 
in Fig. 4, are plotted for the moderate initial rates. 

In Fig. 5, we analyze the normalized outage factor and evaluate the theoretical results of 
Theorem 2. Considering quasi-static conditions, M = 1, At = 1 and f = 5dB, Fig. 5a shows 
the outage probability versus the product of the number of transmit and receive antennas. Also, 
Fig. 5b demonstrates the normalized outage factor in Case 1 and compares the results with the 
theoretical derivations of Theorem 2. Finally, Fig. 5c studies the outage probability in Case 2 
and compares the slope of the curves with the normalized outage factor derived in Theorem 2. 
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Here, the results are obtained for the slow- and fast-fading eonditions (T = 2) with i? = 1, M = 
1, = 1 and 0 = 5dB. 

Figure 6 evaluates the effeet of non-ideal PAs and adaptive power alloeation on the performance 
of large MIMO setups. Considering fast-fading conditions with T = 2, Case 2 with large 
(resp. given) number of transmit (resp. receive) antennas and the outage probability constraint 
Pr(Outage) <9.,9 = 10“^, Fig. 6a demonstrates the supported initial transmission rates, i.e., the 
maximum rates for which the outage probability is guaranteed, versus the total consumed power. 
For the non-ideal PA, we set = 30 dB,'i9 = 0.5, e = 0.65, while the ideal PA corresponds 
to —)-oo,'(9 = 0,e=lin (l25l) . The figure demonstrates the simulation results while, with 

the parameter settings of the figure, the same (with high accuracy) results are obtained if the 
supported initial rates are derived analytically according to (l26l) (Also, see Fig. 2 for the tightness 
of approximations in Case 2). Finally, assuming slow-fading conditions and Case 2 with A, = 1, 
9 = 10“^, M = 2, Fig. 6b compares the required number of transmit antennas in the scenarios 
with and without adaptive power allocation between the HARQ retransmissions. Note that, in all 
cases, we have also investigated a wider range of parameters and fading conditions, but because 
the performances of those cases have followed the same trends as the ones shown, we have not 
included those results to avoid unnecessary complexity. According to the results, the following 
conclusions can be drawn: 

• For Cases 1-3 and different fading conditions, the analytical results of Theorem 1 are very 
tight for a broad range of initial transmission rates, outage probability constraints and SNRs 
(Figs. 1-3). Also, in Case 1 (resp. Case 2) the tightness of the approximations increases with 
the number of receive (resp. transmit) antennas (Figs. 1-2). Moreover, the approximation 
scheme of Theorem 1 can accurately determine the required number of antennas in Case 3 
with different values of K. For Case 4 (which is not of practical interest in large MIMO 
setups), we can find the required number of antennas accurately through Theorem 1 when 
(fTTI) is solved numerically via ([8]). As such, the approximations (6) — (c) of (fT5])-(fT^ decrease 
the accuracy, although the curves still follow the same trend as in the simulation results. 
For instance, with different approximation approaches of Case 4, the required number of 
antennas increases with the initial rate linearly, in harmony with the simulation results (Fig. 
4). The tightness of the approximations in Cases 3 (resp. Case 4) increases when the SNR 
decreases (resp. increases). Finally, the scaling laws of Theorem 1 are valid because, as 
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demonstrated in Figs. 1-4, in all eases the analytieal and the simulation results follow the 
same trends (see Theorem 1 and its following diseussions). 

• In all Cases, better approximation is aehieved via Theorem 1 in fast-fading (resp. slow- 
fading) eonditions eompared to slow-fading (resp. quasi-statie) eonditions. This is intuitively 
beeause the eentral limit Theorem provides better approximation in Lemma 1 when the 
number of experieneed fading realizations inereases. 

• The required number of antennas deereases as the outage probability eonstraint is relaxed, 
i.e., 6 inereases, while for different transmission SNRs, there is (almost) a fixed gap between 
the eurves of different outage probability eonstraints (Fig. 1). Also, fewer antennas are 
required when the number of fading realizations experieneed during the HARQ paeket 
transmission inereases. Intuitively, this is beeause more diversity is exploited by the HARQ 
in the fast-fading (resp. slow-fading) eondition eompared to the slow-fading (resp. quasi- 
statie) eonditions and, eonsequently, different outage probability eonstraints are satisfied 
with fewer antennas in the fast-fading (resp. slow-fading) eonditions (Fig. 2). However, 
the gap between the system performanee in different fading eonditions deereases with the 
number of antennas (Fig. 2). 

• In Cases 2 with many antennas only at the transmitter, the performanee improvement 
becomes limited when the number of transmit antennas reaches, say. A) > 70. In the 
meantime, considerable improvement is achieved by adding more antennas at the receiver 
(Fig. 2. Also, similar arguments hold for Case 1 although not demonstrated in the figures). 

• The HARQ reduces the required number of antennas significantly (Fig. 3). For instance, 
consider the quasi-statie conditions, the outage probability constraint Pr(Outage) < 10“^, 
A) = 5, 0 = 5 dB and the code rate 20 npcu. Then, the implementation of HARQ with a 
maximum of M = 2 retransmissions reduces the required number of receive antennas from 
95 without HARQ to 15 (Fig. 3). Moreover, the effect of HARQ increases with the number 
of transmit/receive antennas (Fig. 3). 

• Different outage probability requirements are satisfied with relatively few antennas. For 
instance, consider a SIMO setup in quasi-statie conditions and M = 1,0 = 5dB. Then, 
with an initial rate R = 3 npcu, the outage probabilities Pr(Outage) < 10“^, 10“"^ and 10“® 
are satisfied with 16,18 and 20 receive antennas, respectively (Fig. 5a). These numbers 
increase to 31, 35, and 38 for = 4 npcu (Fig. 5a). 
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• The normalized outage factor, i.e., the negative of the slope of the outage probability curve 
versus the product of the number of antennas as the number of antennas increases, follows 
the theoretical results of Theorem 2 with high accuracy (Figs. 5b and 5c). Moreover, 
the number of fading realizations experienced during the packet transmission increases 
the normalized outage factor linearly (Fig. 5c. Also, see Theorem 2 and its following 
discussions). 

• The inefficiency of the PAs affects the performance of large MIMO setups remarkably. For 

instance, with the parameter settings of Fig. 6a and i? = 10 npcu, = 2, the inefficiency 
of the PAs increases the consumed power by ~ 11 dB (Fig. 6a). However, the effect of the 
PAs inefficiency decreases with the SNR which is intuitively because the effective efficiency 
of the PAs = e(^^)’^ is improved at high SNRs. On the other hand, adaptive power 

allocation between the HARQ retransmissions reduces the required number of antennas 
marginally (Fig. 6b). Therefore, considering Fig. 6b and the implementation complexity of 
adaptive power allocation, non-adaptive power allocation is a good choice for large MIMO 
systems. 

B. On the Effect of Spatial Correlation 

Throughout the paper, we considered IID fading conditions motivated by the fact that the 
millimeter-wave communication, which will definitely be a part in the next generation of wireless 
networks, makes it possible to assemble many antennas close together with negligible spatial 
correlations m, a. However, it is still interesting to analyze the effect of the antennas spatial 
correlation on the system performance. For this reason, considering Case 2 with = 1, Fig. 
7 demonstrates the required number of antennas in the spatially-correlated conditions where, 
denoting the transpose operator by ()^, the successive elements of the channel vector H = 
[hi,..., follow 

hi = fShi-i -f y/l — ffw, w ~ CMifS, 1), /io ~ CM{f, 1). (30) 

Here, ^ is a correlation coefficient where /) = 0 (resp. /3 = 1) corresponds to the uncorrelated 
(resp. fully correlated) conditions. This is a well-established model considered in the literature 
for different applications, e.g., Il48ll . 
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As shown in the figure, the effeet of the antennas spatial eorrelation on the required number 
of antennas is negligible for eorrelation eoeffieients of, say, jS < 0.4. This is in harmony with, 
e-g-, [SHl, IfSOll whieh, with different problem formulations/metries, derive the same eonelusion 
about the effeet of the antennas eorrelation on the system performanee. Then, the sensitivity to 
the spatial eorrelation inereases for large values of the eorrelation eoeffieients, and the required 
number of antennas inereases with /3. However, the important point is that the eurves follow the 
same trend, for a large range of eorrelation eoeffieients (Fig. 7). Thus, with high aeeuraey, the 
same sealing laws as in the IID seenario also hold for the eorrelated eonditions, as long as the 
eorrelation eoeffieient is not impraotieally high. Moreover, we observe the same eonelusions in 
the other eases, although not demonstrated in the figure. Finally, it is worth noting that, as shown 
in m, for moderate/large number of transmit and/or reeeive antennas and with appropriate mean 
and varianee seleetion, the aeeumulated mutual information of the eorrelated MIMO setups 
follows Gaussian distributions with high aeeuraey. Therefore, one ean use (ISTll and the same 
proeedure as in our paper to derive elosed-form expressions for the required number of antennas 
in the spatially-eorrelated MIMO-HARQ systems. 

VI. Conclusion 

This paper studied the required number of antennas satisfying different outage probability 
eonstraints in large but finite MIMO setups. We showed that different quality-of-serviee require¬ 
ments ean be satisfied with relatively few transmit/reeeiver antennas. Also, we derived elosed- 
form expressions for the normalized outage faetor whieh is defined as the negative of the slope of 
the outage probability eurve plotted versus the produet of number of antennas. As demonstrated, 
the required number of antennas deereases by the implementation of HARQ remarkably. The 
effeet of the antennas spatial eorrelation on the required number of antennas is negligible for 
small/moderate eorrelation eoeffieients, while its effeet inereases in highly eorrelated eonditions. 
Finally, with the problem formulation of the paper, the performanee of the large MIMO systems 
is sensitive (resp. almost insensitive) to the power amplifiers ineffieieney (resp. adaptive power 
alloeation between the HARQ retransmissions). Performanee analysis in the eases with other 
HARQ protoeols is possibly an interesting extension of the work presented in this paper. 
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Figure 1. The required number of receive antennas vs the initial transmission rate R (Case 1: large N,, given N,). Outage 
probability constraint Pr(Outage) < 6 (6 = 10“"^ or 10“^), slow-fading conditions, M = 2, and Nt = 1. 
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Figure 6. Subplot (a): On the effect of non-ideal PAs. Supported initial transmission rate vs the consumed power 10 log 0 ( 0 “"*). 
Case 2: (large Nt, given Nr), outage probability constraint Pr(Outage) <0,6 = 10““^, T = 2, M = 2, fast-fading conditions. 
In the cases with non-ideal PAs, we set e = 0.65, t? = 0.5,0““ = 30 dB. Subplot (b): On the effect of adaptive power 
allocation. The required number of transmit antennas vs the initial transmission rate R (npcu). Case 2: (large N, given Nr), 
outage probability constraint Pr(Outage) < 6 with 6 = 10“®, Nr = \,4> = —5 or 0 dB, M — 2, and slow-fading conditions. 
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Figure 7. The required number of antennas in different spatially-correlated conditions. Case 2: (large Nt, given Nr), outage 
probability constraint Pr(Outage) < 6 with 6 = 10“"^, M — 1, quasi-static conditions, and Nr = 1. 
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